1 

Multilevel Coding over Two-Hop Single-User 

Networks 



Vahid Pourahmadi, Alireza Bayesteh, and Amir K. Khandani 

Coding & Signal Transmission Laboratory (www.cst.uwaterloo.ca) 

Dept. of Elec. and Comp. Eng., University of Waterloo 

. Waterloo, ON, Canada, N2L 3G1 

Q . Tel: 519-725-7338, Fax: 519-888-4338 

O 

' e-mail: {vpourahm, alireza, khandani} @cst.uwaterloo.ca 



Abstract 

In this paper, a two-hop network in which information is transmitted from a source via a relay to a destination 
^ ■ is considered. It is assumed that the channels are static fading with additive white Gaussian noise. All nodes 

are equipped with a single antenna and the Channel State Information (CSI) of each hop is not available at the 
^ corresponding transmitter. The relay is assumed to be simple, i.e., not capable of data buffering over multiple 

, coding blocks, water-filling over time, or rescheduling. A commonly used design criterion in such configurations 

■ is the maximization of the average received rate at the destination. We show that using a continuum of multilevel 

^ : ... . ... „ . ... . _ ™. »a ... ... 

' optimum in this setup. In addition, we present a scheme to optimally allocate the available source and relay powers 

, to different levels of their corresponding codes. The performance of this scheme is evaluated assuming Rayleigh 



X 



fading and compared with the previously known strategies. 

I. INTRODUCTION 

In recent years, relay-assisted transmission has gained significant attention as a powerful technique to 
enhance the performance of wireless networks. The main idea is to employ some extra nodes (relay nodes) 
in the network to facilitate the communication between the terminal nodes. The concept of relaying was 
first introduced by Van der Meulen in [1] and is defined as a scheme to improve the coverage/reliability 
of a wireless network. For instance, relays are usually deployed in networks when the direct link between 
the source and the destination is either blocked or has a very poor quality. The term two-hop network 
usually refers to such a network configuration in which there is no direct link between the source and 
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the destination, and one relay node assists the transmission of data between the end terminals, see Fig. [U 
Two-hop networks have been implemented widely in different applications, including TV broadcasting 
and satellite communications. 

Following the introduction of relay channel in [1], Cover and El Gamal introduce two different coding 
strategies for single relay networks [2]. In the first strategy, known as "Decode and Forward" (DF), the 
relay decodes the transmitted message and cooperates with the source to send the message in the next 
block. Instead of decoding, in the second strategy, the relay compresses the received signal and forwards it 
to the destination in the next block. The terms "Compress and Froward" (CF) or "Quantize and Forward" 
(QF) usually refer to this transmission scheme. Besides DF and CF, in some recent results, [3]-[6], the 
authors investigate another transmission scheme called "Amplify and Forward" (AF) for the Gaussian 
relay network. In this strategy, without decoding the information, the relay amplifies the received signal 
and retransmits it to the destination. 

Knowing these schemes, the performance of relaying is analyzed for different network topologies. For 
instance, considering a single-relay network, authors of [5] and [6] derive a single-letter expression for 
the maximum achievable rate of AF relaying using a simple linear scheme (assuming frequency division 
and AWGN channel). As another example, [4] shows that AF relaying achieves the network capacity 
of Gaussian parallel single-antenna relay network. The extension of [4] to the case of multiple- antenna 
Rayleigh fading networks is presented in [7] and [8]. The first capacity result for relay networks is obtained 
in [2], where the authors prove the optimality of DF strategy in a single-relay network when the received 
signal at the destination is a degraded version of the relay received signal. Clearly, the degradedness 
condition holds in the two-hop setting. Thus, DF would be the optimal relaying scheme for two-hop 
networks. Indeed, most of the results in the literature on relay networks either assume static channels 
between the nodes or perfect knowledge of the Channel State Information (CSI) at both end nodes of 
each link, for the case of fading channels. 

Recently, some papers discuss different transmission schemes over relay networks when the CSI is 
not available at the transmitting nodes, where most of them focus on Diversity-Multiplexing Trade-off 
(DMT) [9]-[14]. Obviously, for such settings, the ergodic capacity is not defined, however, the outage 
capacity is defined as the maximum rate decodable with a given probability [15]. From the throughput 
maximization point of view, the goal is to propose a scheme to maximize the average data rate received 
at the destination. The simplest form of such a problem is to find the optimal transmission strategy for a 
one-hop single-user network when CSI is available only at the receiver and the channel has quasi-static 
fading characteristic. In a pioneering work, Shamai has addressed this problem by substituting the receiver 
by a continuum of virtual receivers, each corresponding to a specific realization of the channel gain [16]. 
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Relying on the resulting degraded broadcast channel, [16] shows that an infinite-level coding scheme with 
a proper power allocation among the different levels of code maximizes the destination's average data 
rate. 

There are several extension for this work; in [17]-[19] the authors try to find the optimal transmission 
strategy when partial information is available at the source. [20] suggests the application of multilevel 
coding in a multicast network with some QoS constraints and derives the optimum throughput-coverage 
trade-off for such a network. [21] combines multilevel coding scheme with Hybrid Automatic Retrans- 
mission Request (HARQ) and shows that this approach results in high throughput and low latency in a 
point-to-point link. The Multiple Input Multiple Output (MIMO) extension of [16] is also discussed in 
[22], [23], and [24]. 

Considering the two-hop network with no CSI at the transmitter side of each link, Steiner et al in 
[25] propose different transmission schemes, assuming that the relay node has a power constraint and 
is simple, i.e., is not capable of data buffering over multiple coding blocks or rescheduling. To find the 
optimal transmission scheme, they study different multilevel coding schemes when the relay operates in 
different modes, including DF, AF, and CF. As discussed in [25], due to the high complexity of the 
infinite-level DF coding scheme, they have only considered a finite level code in their proposed DF 
strategies. Comparing the results of [25], it turn out that AF strategy outperforms all other investigated 
transmission schemes, specifically in the high SNR regime. However, as concluded in [25], although AF 
has the best performance among the other strategies, the optimality of AF scheme is not implied. In 
fact, since the general infinite-level DF coding scheme remains unsolved, there is always a question of 
whether or not DF relaying achieves a higher expected rate as compared to AF. The main motivation 
of this work is to answer this question. 

In this work, we study the performance of DF relaying scheme which uses infinite-level codes at the 
source and the relay in the same setup as in [25]. To this end, we first prove that the infinite-level code 
and DF relaying is indeed the optimal two-hop transmission strategy for maximizing the destination's 
average data rate. We further propose an algorithm to determine the optimum power allocation for each 
indefinite-level code. Numerical results are presented to verify that the proposed scheme outperforms the 
AF strategy discussed in [25]. 

The organization of this paper is as follows: First, section HI] describes the two-hop system model. A 
brief review of some previous results on single-hop links are presented in section Hill Formulating the 
DF multilevel coding scheme in subsection IIV-A[ we prove the optimality of this scheme in subsection 
IIV-BI Afterwards, in section llV-Ci we present a procedure to actually determine the optimal infinite-level 
code parameters at the source and the relay. Numerical results and comparison with other schemes for the 
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A K) 

Source Relay Destination 

Fig. 1. Two-hop Network Model 

special case where both links are Rayleigh fading is presented in section |Vl Finally, section |Vl] concludes 
the paper. 

Throughout this paper, we represent the expectation operation by E[.]. The notation log(.) is used for 
the natural logarithm and rates are expressed in nats. We denote fy{.) and Fy(.) as the probability density 
function (pdf) and the cumulative distribution function (CDF) of random variable y. 

II. System Model 

In this paper, we investigate the performance of a two-hop network, in which a relay assists the 
transmission of data between a source and a destination. As Fig. \T\ shows, in a two-hop network, the 
destination can solely receive data via the relay. All nodes are assumed to have single antennas. It is 
assumed that the source has no information about either of the channel gains, the relay knows only 
the channel gain between itself and the source, and the destination knows both channel gains. This CSI 
assumptions are indeed practical, since the receiver of each hop can evaluate its immediate channel gain by 
measuring the pilot signal sent from the corresponding transmitter. In addition, the receiver can measure 
the equivalent channel (the source to the destination) if the relay forwards the pilot signal of the source 
towards the destination. Having the equivalent channel gain and the relay to destination channel gain, the 
destination can find the source to relay channel gain as well. 

It is also assumed that the source and relay both know the fading power distribution of both links. By 
/r(7) and friv) we denote the probability density function of fading power in the first and second hop, 
respectively. The channels of both links are assumed to be quasi-static, i.e., they are chosen randomly 
(based on their corresponding pdf) at the start of the transmission and remain fixed for the whole 
transmission. 

Source and relay are assumed to have power constraints Pg and P^, respectively. Furthermore, they are 
assumed to be capable of constructing infinite-level codes. However, the relay is assumed to be simple, 
i.e., is not capable of buffering any unsent data, rescheduling the untransmitted data from the previous 
blocks. It is also assumed that the relay can not perform water- filling of its power over multiple blocks 
and has a maximum power constraint in each block. Therefore, it can only retransmit the data which 
it has just received from the first hop. Successive decoding is used as the decoding procedure at both 
destination and relay (in case the relay wants to decode the data). 
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III. Multilevel Coding Scheme for Single-hop Networks 



In this section, we review some previous studies on the application of the multilevel codes for a single- 
hop network in two scenarios: 1) No CSI is available at the transmitter, which is known as "broadcasting 
approach" [16] and 2) no CSI is available at the transmitter and in addition, the transmission data rate 
can not exceed a certain value [25]. 

A. Single-Hop Broadcasting Strategy 

The optimum scheme for a single-hop link, in case that the transmitter knows the fading power for 
each block is to design a single-level code with the rate of log(l + IP) for that block. Here, P denotes 
the normalized transmission power, i.e, the equivalent transmission power if the noise power is equal to 
one and / is the fading power for that specific block. Therefore, the average achievable rate for this setup 
would be Rerg = El [log(l + IP)] [26]. 

Although designing a single-level code is optimum in the above scenario, it can not be applied for a 
transmitter which does not have access to the channel state information. For such a scenario, Shamai has 
introduced a technique called broadcasting strategy [16]. In this technique, the transmitter sends the data 
through infinite levels of a superposition code. Then, conditioned on the channel state, i.e., the fading 
power, the receiver decodes up to a certain level of the code. Therefore, the total receiving rate for each 
channel realization, say I, can be evaluated as: 





where dR{a) represents the differential rate transmitted over level 'a' of the code. Defining p{l)dl as the 
power assigned for the level, dR{l) is given by: 



where /(/) = p{a)da and the second line follows from the assumption of infinitesimal rate assignment 
which is shown to be optimal in [20] . The aim is to find a p(.) such that the average data rate at the 
destination is maximized while the source power constraint (P) holds. It means that the power should be 
assigned to different code levels such that: 




(1) 




lp{l)dl 



(2) 



! + //(/)' 



oo 



max 




(3) 







oo 



s.t. 
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where /(.) shows the pdf of the channel fading power. 

This maximization problem has been studied and solved using calculus of variations technique (see 
[16] for details of the proof). Here, we only mention the final solution as: 



P 1<Iq 
li<l 



where F{1) = f{a)da is the CDF of the fading power, lo and h are determined such that they satisfy 
/*(/o) = P and = 0, respectively. Clearly, the optimum power assignment can be determined by 

p*{l) = The maximum destination's average data rate will be: 

Rl, = / dlf{l)R*{l), (5) 



where R*{1) is obtained by setting p{l) = p*[l) in ([U and The total rate of the final superimposed 
code can be evaluated by: 

poo 

R*F= dR*{l). (6) 
Jo 

Note that this value only depends on the fading power distribution and the transmitter power. 

Finally, it is important to mention that the above scheme is indeed the optimal transmission scheme 
for a fading link when the source does not have the CSI. It is due to the result of [27] which proves that 
the multilevel coding maximizes any weighted sum-rate of a degraded broadcast channel. Therefore, 
considering the equivalent broadcast model for the single-user fading channel, the multilevel coding 
achieves the optimal destination's average data rate. 

B. Single-Hop Rate-Limited Broadcast Strategy 

An interesting extension of the above broadcast strategy is designing a multilevel code for a source 
with available data rate limited to i?j„, i.e., the transmission rate should be less than [25]. 

The problem formulation is similar to subsection IIII-Ai except it has one more constraint on the 
transmission rate. The modified optimization problem will be, [25]: 

oo 

max Ray = [ dlf{l)R{l) (7) 
p(-) J 



oo oo 



s.t. j p{l)dl = P, and J dR{l) < R^ 





where the second condition ensures that the transmission rate remains less than the available source data 
rate. 
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Reference [25] uses constrained calculus of variations to solve this problem. More precisely, initially 
the authors in [25] substitute the first constraint by two end-point conditions of J(0) = P and /(oo) = 0. 
Then, using p*{l) = — '^^J''^ and Lagrangian multiplier, (|7]) is reconstituted as a variations problem (See 
[25] for more details). It turns out that the optimum distribution function is as follows: 

^ P l<lo 



r(0 



7(ljp h < I' < I'l 

h<l 



(8) 



where F{1) = J f{a)da. Iq and li are determined as a function of A and satisfy /*(/o) = P and = 0, 

— oo 

respectively. Finally, A is computed such that the transmission rate constraint holds. 



As an example, for the case of Rayleigh fading channel, i.e., F(l) 
is as follows [25]: 

P I <lo 

r{i) = 



the optimum distribution 



A 





I lo<l <li 
h < I 



h = l-WU-Xe) 



(9) 



(10) 



where Wl{x) is the Lambert W-function and finds lu such that we"^ = x. The values for A and Iq are also 
determined by solving the following system of equations: 

Rin = 21og(/i)-/i-(21og(/o)-/o) ^^^^ 
/*(/o) = P 

One important observation is that in the case of Rm > R*p {R*p is defined in equation Q) the above 
rate-limited problem will be simplified to the original problem in subsection IIII-Ai Although this statement 
can be verified mathematically, its intuitive explanation would be insightful. It is obvious that if i?j„ tends 
to infinity, the rate condition is always satisfied (will not be an active constraint). Therefore, maximization 
problem of © relaxes to the form of equation ([3]). Moreover, IIII-AI shows that the source, in the optimal 
transmission of the original problem (without rate constraint), feeds the channel with a rate equal to Rp. 
Hence, it turns out that even though the available rate at the source is more than R*p, the source needs 
to transmit only R*p bits of information in each block. Therefore, if Rin > R*p, the solutions of the two 
optimization problems of ^ and ^ are equal. 

At the end, note that the resulted multilevel code is also the optimal transmission scheme for the 
rate-limited one-hop set-up. To prove, we should use the broadcast equivalent structure of the network. 
Indeed, the capacity region of the rate-limited degraded broadcast network is the intersection of the 
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Single-Hop Rate Limited Capacity Region 
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Fig. 2. Single-Hop Rate Limited Capacity Region 

original degraded broadcast channel capacity region and the region below the surface associated with the 
transmission rate constraint. For illustration, the solid line in Fig. [2] shows a typical capacity region of a 
rate-limited two-user degraded broadcast channel. In this figure, dotted and dashed lines depict the capacity 
region of the original degraded broadcast channel (without rate limitation) and the line representing the rate 
constraint, respectively. The destination's average data rate can be evaluated as the sum of the received 
data rate for each channel state / (-R(O) times the probability of occurring that specific channel state 
(/(/)). Clearly, this value is maximized on the boundary of the resulted capacity region (the solid line). 
Therefore, the optimal point would be either over the capacity region without rate limitation (Arc AB) 
or the end point of the rate limitation line (Point C). Since AB is a part of the original capacity region, 
multilevel coding is the optimal scheme to achieve any point on AB. Furthermore, point C is archived 
using a single level code which is again a special case of multilevel codes. This proves the optimality of 
multilevel coding scheme for rate limited scenarios. 

IV. Multilevel Coding Scheme for Two-hop Networks 

As an extension of the one-hop set-up, [25] addresses the problem of maximizing the destination's 
average data rate in a two-hop network, where there is no direct link between the source and the destination. 
In reference [25], several schemes have been studied, including broadcasting strategy with AF relaying, 
and DF relaying with finite level broadcasting at the source and the relay. Infinite level codes with DF 
relaying is also addressed in [25]. However, the performance of this method remains as an open problem 
there. 

In this section, we will first describe the infinite level DF strategy in details. Then, in subsection IIV-B[ 
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we prove the optimality of this scheme. Finally, subsection IIV-CI presents an algorithm to optimally design 
such an infinite level code. 



A. DF Infinite-Level Codes for Two-hop Networks 

Based on the system model, each transmission block of the infinite level DF strategy consists of the 
following two steps: 

1) In the first phase, the source allocates its power among different code levels with the power 
distribution function Psi-)- Of course, ps{.) should satisfy the power constraint ps(a)da = Pg. 
Then, based on the source-relay channel fading power, say 7, the relay is able to decode up to the 
level 7 of the transmitted data. Thus, the relay received rate is: 



where /s(a) = ps{a)da. 

2) In the second phase, the relay should transmit the data to the destination. As noted earlier, in this 
work, we only focus on simple relays which can neither buffer any of the previously received data 
nor do any scheduling tasks. As a results, these relays have two features which seem obvious but 
have important effects on the code design. To illustrate, consider a case in which the relay has 
decoded Rr{'j) bits of the transmitted data. It turns out that, firstly, the relay can not transmit with 
the rate greater than Rri'j). Secondly, if the relay transmits with the rate R2, R2 < Rr{l), the rest 
of the data (Rri'j) — R2) can not be stored and should be discarded. Consequently, the relay, in each 
transmission block, should choose the optimal power distribution of the multilevel code such that 
it satisfies the relay total power constraint (Pr). Meanwhile, the relay should keep the transmission 
rate below its received data rate in that block (Rr{j)). 

Defining pr{.\Rr{'j)) as the power distribution of each code level at the relay conditioned on the 
input rate of Rri'j), we can summarize these conditions as: 

a) Power constraint at the relay: Vi?r(7) : 

00 

/ Pr{a\Rr{'y))da = Pr. 


b) Available rate constraint at the relay: Vi?r(7) : 

00 

/ nOT^yy ^ ^^■(^)' ^here i?.(7) is defined by m. 



Clearly, the relay requires to know pr{a\Rr{'y)) for all possible values of Rri'j). 

Transmitting a multilevel code on the relay-destination link, the destination is able to decode up 

to a certain level 'f'. Here, 'f ' denotes the fading power of the second link. Therefore, for each 
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Rrij), the received data rate at the destination can be written as: 

RMM-,)) = /log (i + ;:"':'fi'"y; ) - / ttttw^^"- 

J V I + alrialRri-f)) J J 1 + a/r.(a|i?r(7)) 



Indeed, for successful decoding of the signal, the destination should know the power allocation 
strategy of the relay. This information can be obtained through the knowledge of the source to the 
relay channel gain. 

Given these, we are now able to formulate the two-hop optimization problem. Similar to the single-hop 
scenario, we want to maximize the average data rate received at the destination. Assuming /r(7) and 
fri"^) as the probability density functions of the fading power in the source-relay and relay-destination 
links, respectively, the destination's average data rate can be written as: 

E[Rd] = Er{Er[RAv\Rr{im (14) 



Therefore, we obtain the final optimization problem as follows: 





max / I fr{i)h{v)RMRr{l))dv d-i (15) 

Ps(-), Pr{.) 







s.t. J ps{a)da = Ps, 



oo 

Vi?r(7) : J Pr{a\Rr{l))da = P, 





oo 



apr{a\ 


\Rr{'j))da 


1 + alr{a\ 


\Rril)) 



Note that the above optimization problem is similar to the one derived in [25]. However, in [25], the last 
constraint (relay rate limitation) is stated as an equality. In fact, it is more accurate to formulate the rate 
limitation by an inequality constraint instead of equality. It is due to the fact that the relay may not have 
to send all information it receives from the first hop to achieve the optimal performance. In other words, 
the last constraint of equation (fT5l) lets the relay to discard some of its received data if it wants to do so. 
For instance, this may happen when the relay receives data rate higher than its corresponding R*p (R*p is 
defined in equation In such a scenario, the relay only uses R*p bits of the received information and 
ignores the rest. 
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Fig. 3. Virtual Two-Hop Broadcast Channel Model, k = 3 

B. Optimality of Two-Hop DF Multilevel Coding 

The main focus of this section is to show that the multilevel coding approach combined with decode 
and forward (DF) relaying maximizes the average data rate at the destination of a two-hop network. 

To start, let us emphasize that according to the two-hop structure of the network, all information received 
by the destination should be first passed through the relay and there is not any direct link between the 
source and the destination. As a result, the destination received signal is always a degraded version of 
what has been received at the relay. In other words, no information can be decoded by the destination 
unless it has been decodable at the relay. Thus, it can be concluded that decode and forward is the optimal 
relaying scheme for two-hop settings. 

Knowing the optimality of DF, similar to the single-hop network (See [22]), we model both first and 
second fading hops by infinite number of virtual relays and virtual users, respectively. More precisely, we 
substitute the relay with infinite relays, each has a constant channel gain which corresponds to a specific 
realization of the first hop channel. These virtual relays constitute a degraded set. If we assume that the 
channel fading power is selected from a set of discrete values {ai, 02, ■ " " ) O'k^ there would be k virtual 
relays in the network. Without loss of generality, we assume ai < 02 < ■ ■ • < a^. Of course, this model 
is accurate only if k tends to infinity. In a similar way, the second hop can be modeled with k virtual 
users for each of the virtual relays; thus, in total there would be virtual users in the network. Uij 
denotes the virtual user which is associated to the virtual relay Ri, and its channel gain is a^. Moreover, 
Bij represents the decodable rate at this node. For illustration purpose. Fig. |3] depicts the network model 
for the case that k = 3, i.e., A = {ai, 02, 03}. 

The aim is to find the optimum transmission scheme for the source and the relay. We start from the second 

'Note that here, for simphcity, we assume that the fading levels of both links are selected from the same set. However, the statements of 
optimality holds for the general case. Moreover, this assumption is valid for the case of k ^ 00 
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hop and assume the source uses an arbitrary transmission scheme. Furthermore, let {Bi, B2, - ■ ■ , B^} be 
the rate that each of the virtual relays, i.e., {Ri, i?2, ■ ■ ■ , Rk}, can decode under this transmission scheme. 
Having the data rate Bi at relay Ri, the optimal second-hop strategy is to maximize the destination's 
average data rate for each of the virtual relays. As will be described in subsection IIV-C[ the solution to 
this problem is similar to the result of the subsection IIII-BI and the optimal power distribution, p*Q,{u), 
for each of the virtual relays can be determined by dS]), in which = B^. In other words, multilevel 
coding is the optimal strategy for the second hop. 

Note that the successful decoding is only possible if the destination knows the applied p^. (m) (or 
equivalently the value of Bi) at the relay. In our two hop model, this information can be obtained by the 
knowledge of the source to the relay channel gain at the destination. The optimal received data rate for 
each of the virtual users, say Uu, can be evaluated by: 

/up%, (u)du 
, V, , (16) 



where I%^{u) = p*^,{a)da. Let iprU) be probability that the second hop gain is equal to aj. Hence, 
the optimal average rate of the second hop under rate condition Bi would be: 

K 

D*{B,) = Y,i^r{l)R*{B„ T = ai). (17) 
1=1 

Given D*[Bi), the two-hop network average data rate can be written as: 

K 

R = Y,Ml)D*{Bi), (18) 

1=1 

where tpr{l) represents the probability that the first hop channel is in state a^. The goal of the code 
design is to maximize R. As (fTSi) shows, the destination's average data rate is the weighted sum of a 
non-linear function of B^. The domain of acceptable Bi for i E {1,2,...,k} is a convex set which is 
known as the capacity region of the underlying broadcast channel. Moreover, due to the degradedness of 
the virtual relays, all points on the boundary of this capacity region can be achieved by a multilevel code 
[2]. Therefore, to prove the optimality of the multilevel code for the first hop, it is enough to show that 
R is maximized on the boundary of this capacity region. This argument can also be justified if we can 
show that ^ is positive Vz G {1, 2, k}. To show this, we write: 

BR 

^ = M^)D*'{B,) (19) 

where D*(Bi) is defined in equation (flTl) and shows the destination's average data rate when the relay 
has the rate constraint of Bi. By definition, iprii) is a non-negative value. Moreover, from the results 
of subsection IIII-Bl it can be concluded that D*\Bi) is non-negative. In fact, if Bi > R*p, where R*p 
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is defined in equation ©, then D*'{B.i) = and D*'{Bi) > 0, otherwise. Therefore, §§: > for 
Vz G {1,2, k}. Consequently, R is maximized for the values of {Bi}f^-^ which are on the boundary of 
the capacity region of the first hop broadcast channel. In other words, the multilevel coding scheme is the 
optimal transmission strategy for the first hop. This completes the proof for the optimality of multilevel 
coding for the two-hop network. 

C. Optimal Design of Two-Hop DF Multilevel Coding 

Having the optimality of two-hop multilevel coding scheme, in this section, we present a procedure 
in order to solve the two-hop optimization problem introduced in equation (fT5l) . The main difficulty of 
this problem is that, unlike single-hop scenarios (the original and the rate-limited broadcasting cases), 
equation (fT5l) can not be directly solved by variations methods. It is due to the fact that the constraint on 
the second hop rate does not have a fixed value on the right hand side, i.e., it does not have a form of 
isoperimetric problem. For a complete discussion on isoperimetric problem, refer to [28]. To solve this 
problem, we use the following lemma. 

Lemma 4.1: The following maximization problems on two non-negative functions and /2(-)- 

max / / dxdyn{x)}C{fi{.), f2{.),x,y) , (20) 

/i(-)j2(.) Jx=a Jy=c 

and 

pb I'd 

max/ dxn{x) max K:{h{.), f2{.),x,y) dy, (21) 

hi-) Jx=a {f2{-)\ifii-),x)} Jy=c 

where 7i(.) and /C(.) are two known non-negative functions, are equivalent. 

Proof: Let us denote the solution of ^ by (/*(.), /*(.)) and the solution of dZB by (/*(.), /*(.)). 
We can write 

/ dxdyH{x)IC {/;{.), f;{.),x,y) = / dxH{x) [ /C (/*(.),/*(.), y) dy 

x—a J y—c J x—a J y—c 

< / dxH{x) max IC {f^{.), f2{.),x,y) dy 

J a {/2(-)l(/r(-),^)} J.y = C 

< max / dxH{x) max / lC{fi{.), f2{-),x,y) dy. (22) 

/l(-) Ja {/2(-)l(/l(-),^)} Jy = c 

On the other hand, 

max / dxn{x) max ]C{fi{.), f2{,.),x,y) dy ^ / dxn{x) IC {/*{.), f^i-), x,y) dy 

fii-)Ja {f2(-)\ifli-),x)} Jy^c Jx=a Jy=c 

= f f dxdyn{x)IC{fn.),f^{.),x,y) 

J X—a J y—c 

- Vl^]^,, / dxdyH{x)lC{fi{.),,f2{-),x,y). 

/l(-)./2(-) J x=a J y=c 

(23) 
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Combining (|22l) and (|23l) . lemma is proved. 

Using this lemma and noting that /r(7) > 0,V7 , we can reform E[Rd] in (fT5] ) as follows: 



oo oo V 



max E\Rd\ = max /" d-y fril) max /" friv) [ - — ^^^^ \ da, 







max /"rf7/r(7) max [ dv fr{v) [ -^^^^^^:f^^da, (24) 



where the outer maximization is subject to: 

oo 

j ps{a)da = Ps, (25) 



and the constraints of the inner problem are as follows: 

oo 

yRrij) : / pria\Rr{-f))da = Pr, (26) 





oo 



Vi?.(7):/TfS|^f^rfa<i?.(7). (27) 



In (|24|) . (a) follows from the fact that i?r(7) can be determined with the knowledge of 7 and ps{.) and 

apr{a\Rr{'r)) 
l+a/r(a|Rr(7)) 



the dependence of the term / dv friv) J t+^/7^rt\)'^'^ Ps(-)i7 through i?r(7)- 





Given (|24|) -(|27]). in the following two subsections, we will discuss how this two-step maximization 
problem can be solved using Euler's equations [29]. 

1) Relay-Destination Link Optimization Problem 

Receiving i?r(7) bits from the first hop, the aim of the relay is to maximize the average data rate received at 
the destination. In fact, if the input rate changes, the relay should modify its power distribution, accordingly. 
However, the knowledge of the input rate (Rr{'y)), the relay total power, and the pdf of the second hop 
fading power is sufficient for determining the optimum power distribution function, p*(.|i?r(7))- It is 
evident that the optimum power distribution function, can be completely determined by evaluating 

p*.{.\Rr{j)) for all values of Rr{'j). pl{.\Rr{'j)), itself, is the solution of the following problem: 





00 



apr{a\ 


\Rr{'j))da 


1 + alr{a\ 


\Rrh)) 



h{Rrh)) = max / dv Mv) / ^^"^ ' . 7. (28) 

Pr(.|i?r(7)) , < „ „ 



.t. J pr{a\Rr{''y))da = Pr, 



00 

r aMa|i?.(7)) 
J 1 + aIria\Rri-f)) - 
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MaKimum Achievable Average Data Rate at the Destination 
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Fig. 4. Maximum Average Data Rate Received at the Destination 

Note that, in (l28l) . Rr{'j) is a constant; hence, the problem takes the form of the rate-limited broadcast 
strategy problem, subsection IIII-BI Therefore, the optimum solution is: 



i-i.,(zHA-z/.(0 i^^i^i^ , (29) 
li<l 



where Fy{l) = J fr{a)da. Iq and li are determined as a function of A to satisfy /j!(/o) = Pr and = 

— oo 

0, respectively. The optimum multilevel power distribution at the relay can be found by p*(/|i?r(7)) = 
_^£rii^r(7)l pinally, A is computed to satisfy: 

oo 

/ ^ — T I \u I w ^^ = mm{Rr{-r),Rp), (30) 
J 1 + aIr{a\Rr{'y)) 



where R*p is defined by This condition comes from the fact that achieving the maximum average 
rate at the destination requires the relay not to transmit more than R*p bits (refer to the discussion in 
subsection IIII-BI) . 

As an example, we have solved (|28l) for a network in which the second hop can be modeled as a 
Rayleigh fading channel, i.e., Fr{v) = 1 — e~*^. Fig. |4] shows the maximum destination's average data 
rate, h*{Rr{'y)), for different relay input rates (i?,.(7)) and different relay powers (Pr). 
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2) Source-Relay Link Optimization Problem 

Knowing the optimum value for the inner integration, h*{Rr{'y)), (|24l) can be written as follows: 

oo 

E[Rd] = max / d-f fri^)h*iRri^)) (31) 

Psi-) J 



oo 

s.t. J ps{a)da = Pg 



where: 







7 < /, 







p op^^^ /0<7<^1 > 
Jlo l+als{aj ^ ' ' 

J/o l+als(aj ^ I 



(32) 



and Is{l) = J ps{a)da. Iq and li satisfy /^(/o) = Ps and /^(/i) = 0, respectively. As (|32l) suggests, Rr{'j) 
I 

only depends on 7, Ps(.)5 and /s(.)- Remembering ps(0 = '^^'^ write the integrand of (|3TI) as 



G{l,Is,Is) = fr{l)h*{Rr{l, Is, I's))- With this notation, equation (|3TI) takes the form of a fixed end-point 
Calculus of Variations problem and can be solved using Euler's equation, [29], 



c(/,/.,/:) = G,,-^ = o, 



di 



(33) 



where Gi^ = w-^f^, Gi> = and is the derivative of Gr with respect to /. Thus, we have: 



dRr dis 



dRr dl' 



dGdRr_\ fvm'{Rr{lJsJ's)) j j^^da k<l<h 



dRr dig 



otherwise 



(34) 



dGv 
dl 



d^dR^^ I MlW{Rr{l,Is,I's))I 

dRr OR ~ 



l+als{a) 



da Iq < I < li 



otherwise 

mh*'iRr{l, Is, I's)) I T^.da + frm*\Rr{l, Is, 1^)^ 



(0 



+ fril)h*"iRril,Is,rs)) 



l+lls{l) J l+als{a) 




da 



(35) 



lo <l < h 



,(36) 



otherwise 

where h*' {.) and h*" {.) denote the first order and the second order derivative of h*{.), respectively. 
Substituting dMl)-® in the optimal I*{1) is derived. 
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Optimal I (I) for the Rayleigh Fading Two-Hop Nehn/ork 
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Fig. 5. Optimal Is{l) for the Rayleigh Fading Two-Hop Network 



As an example, in the scenario where both source-relay and relay-destination links are modeled with a 
Rayleigh fading channel, i.e., Fr{l) = Fx(/) = 1 — e~', (|33l) can be simplified to: 



h*'{i) 



— a — a Is{a) 



da 



(37) 



-da 



1 + lIsil) J 1 + aIsia) 



where i = Rr{l, h, I's)- To solve (l37l) . we first need to have h*' (i) and h*"{i). Indeed, these functions 
can be numerically evaluated using the results of subsection IIV-C1[ In the next step, we replace Ig 
by /s(2), /s(A^)], corresponding to the amount of interference in each leveo Is{mys are in 

descending order, such that /^(l) = Pg and IsiN) = 0. As a result, we have a nonlinear system of 
equations, i.e., C(m, J^,/^) = 0, m = {1,2, A^} which can be solved numerically. The final solution 
for these N variables shows the optimal interference function, /*(/). As an example. Fig. [5] presents /*(/) 
in the case of Rayleigh fading model for both hops and Pg = Pr = 20dB. Having /*(/), the amount of 
power associated for each code level can be determined by = 



^In fact, we have approximated a continuous variable Is{l) with a discrete A'^-level function, which becomes precise as N tends to infinity. 



18 



V. Numerical Results and Comparison with Other Schemes 

In the previous sections, we have proposed a DF muhilevel coding scheme, which is shown to be 
optimum in the underlying network setup, and derived the optimum source and relay power distribution 
through different levels of code. In this section, we compare the performance of the proposed scheme 
(which is the optimal scheme) with the cut-set upper-bound and two other sub-optimal schemes proposed 
in [25]. 

A. Broadcasting Cutset Bound, Ccutset : 

This bound simply says that the achievable average data rate of a two-hop network can not exceed the 
achievable average rate of any of the single-hop links, i.e., the source-relay and the relay-destination 
links. This is independent of the relay structure and its operation. In other words, the cutset bound 
is an upper-bound on the network throughput when we put no limitation on the relay, i.e., the relay 
is capable of buffering the unsent data or rescheduling the buffered data. Therefore, the gap between 
the performance of the proposed scheme and the cutset upper-bound shows the maximum possible 
gain of having a "complicated" relay instead of a simple one. The cutset bound can be written as: 



Crfii.txp.t — min 



da fr{a)Ri{a), / da fr{a)R2{a) 



(38) 



where Ri{a) (i?2(a)) denotes the rate that the relay (destination) can successfully decode when the 
source (relay) transmits over a channel with fading power equal to 'a'. 
B. Amplify and Forward, AF: 

This is the achievable rate of a two-hop network in which the relay performs the amplify and forward 
(AF) on the received source signal. To design the optimum multilevel power distribution, first, the 
total equivalent channel should be evaluated. In other words, the source-relay and relay-destination 
channels combined with AF relaying can be modeled as one channel with a new probability density 
function. Having this new pdf, the optimum power distribution can be evaluated. Details of the proof 
can be found in [25]. The final result can be written as: 



Raf = I da 

where: 



lo 



f°° , a(l + Prl) 

Fs,{x) = 1 - / dl e"'"^^v^, (40) 

Jr'r 

and = -^Fs^{a). Iq and li axe defined such that Iopt{h) = Ps and loptih) = 0, respectively. 

/opt(/) is presented in [25]. 
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Fig. 6. Destination's Average Data Rate Ps=20 dB, P,.=20 - 32 dB 



C. Outage at the Source, Broadcasting at the Relay, DFi_bs: 

This scheme is another suboptimal strategy that has been studied in [25]. In this case, the source 
uses a one level code, known as the outage approach, and the relay uses the optimal multilevel code. 
The subscript "1 — bs" represents the one-level coding and the broadcast scheme at the source and 
relay, respectively. Clearly, this approach is a special case of the optimum DF broadcast strategy, 
i.e., the proposed scheme. The achievable average rate of this scheme can be computed by: 

RDF,i-bs = max(l-Fr(7.)) /cia(l-FT(a))-^^^ (41) 
7s.Pr(.) J i + alr{a) 



oo 

f UPr(u)du , , 

/ 1 V ^ =log(l + P.7s). 

J I + Ulr{u) 


Figures ^ and ([7]) represent the destination's average data rate at the destination versus the relay power 
Pr for the proposed scheme, as well as the AF and DFi_bs schemes, where Pg = 20dB and P, = 30dB, 
respectively. The upper-bound Ccutset is also depicted in both figures. 

As expected, the proposed DF strategy (the optimal scheme) outperforms the AF and DFi_bs schemes. 
Note that, the superiority of the proposed scheme over DFi_bs is obvious since DFi„bs is a special case 
of the proposed scheme. The important observation in these figures is that the infinite multilevel DF 
strategy is strictly superior to the AF strategy, which was previously the best known scheme for this 
setup at the high SNR [25]. However, as the SNR at the relay side Pr increases, the performance of AF 
approaches the optimal performance. Furthermore, as Pr increases, the proposed scheme approaches the 
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Fig. 7. Destination's Average Data Rate Ps=30 dB, Pr=30 - 42 dB 



cutset bound which means that for high values of Pr the relay does not need to be "complicated". Another 
observation from these figures is that, as decreases, DFi_bs approaches the optimal performance. This 
can be explained as follows: when the power of the relay is much smaller than the source power, the 
relay-destination link limits the performance. Therefore, even using a one-level code at the source is 
sufficient to deliver an average rate of R*p to the relay, which is the maximum rate that relay can transmit 
to the destinationy. 

VI. Conclusion 

In this paper, a two-hop network in which the data is transmitted from the source node via a single 
relay to a destination node was considered. It was assumed that the knowledge of the channel for each 
hop is not available at the corresponding transmitter. The relay was assumed to be simple, i.e., not capable 
of data buffering over multiple coding blocks, water-filling over time, or rescheduling. For this network 
setup, we proposed an infinite-level coding scheme at the source and the relay. It is shown that this scheme 
in conjunction with the Decode and Forward (DF) relaying is indeed the optimal strategy for maximizing 
the average data rate received at the destination. We also proposed an algorithm to find the optimum 
amount of power which should be assigned to each code level at the source and relay. The optimality of 
the DF multilevel coding strategy is also verified through numerical results by showing its superiority 
over the Amplify and Forward (AF) scheme, which was previously the best known scheme for the high 
SNR regime. 

^Note that even if relay receives a rate more than Rp, the extra rate should be discarded. 
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